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Abstract 

The first step in the counting operator analysis of the spectrum of 
any model Hamiltonian H is the choice of a Hermitean operator M 
in such a way that the third commutator with H is proportional to 
the first commutator. Next one calculates operators R and E) which 
share some of the properties of creation and annihilation operators, 
and such that M becomes a counting operator. The spectrum of H is 
then decomposed into multiplets, not determined by the symmetries 
of H, but by those of a reference Hamiltonian H re f, which is defined 
by H IC { = H — R — W , and which commutes with M. Finally, we 
introduce the notion of stable eigenstates. It is shown that under 
rather weak conditions one stable eigenstate can be used to construct 
another one. 

Keywords Hubbard model, Jaynes-Cummings model, transverse Ising model, su- 
persymmetry, multiplets, level crossing. 
PACS 71.10.Fd 

In the literature many attempts are found to generalise the notion of creation 
and annihilation operators. Some of these were introduced in the context of Bo- 
goliubov's notion of quasi-particles — see for instance pQ. Others are related to 
the method developed by Darboux in the nineteenth century to find new solutions 
of non- linear equations (see for instance [2]). Lowering and raising operators [3] 
determine recurrence relations and generate a Lie algebra (see for instance [1]). 
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The operators R and w introduced below share some of the properties of these 
operators. In the case of the harmonic oscillator they coincide with the usual 
annihilation and creation operators, up to some additive constants. But in other 
models they can be used to guarantee the existence of a second eigenstate, with 
a different eigenvalue. When applying the technique to a 4-site Hubbard ring 
OE], the spectrum of eigenvalues decomposes into multiplets, some of which have 
a fermionic appearance, in contrast with the bosonic spectrum of the harmonic 
oscillator. 

The one-dimensional Hubbard model has become very famous by the work of 
Lieb and Wu [TJ [8] who used the Bethe ansatz to study the eigenvalues in the 
thermodynamic limit. On the other hand, the energy levels of the 6-site Hubbard 
ring could be determined in a reliable way [9] by making use of the symmetries of 
the Hamiltonian. In the latter work, energy level crossings were observed which 
seemingly were not in agreement with the non-crossing rule of Wigner and von 
Neumann [10] • Their origin was explained in The multiplets considered here 
and in [6] are not determined by the symmetries of H but by those of the reference 
Hamiltonian H Ic f = H — R — Fv . In particular, we do not require integrability 
of the model Hamiltonian H. For instance, our results are also applicable to the 
Hubbard model in dimensions larger than one. 

Consider a quantum Hamiltonian H together with a Hermitean operator M ^ 
0, not commuting with H, such that 

[[[H,M],M],M] = 7 2 [H,M], (1) 

with real 7^0 [12]. Note that such an operator M does always exist. Indeed, let 
M = E where E is any orthogonal projection operator (i.e. E = E' = E 2 ). Then 
(H]) is satisfied with 7 = 1, as is readily verified. The present work generalises this 
observation. 

The relation ([1]) is trivially satisfied if 

[[H,M],M] = 1 2 H. (2) 

This is the case of the Pauli spin. The Hamiltonian is H = ^Hujo<j z . A short 
calculation gives [[H, a x ], a x ] = 4H. Hence, (JTJ) is satisfied with M = a x and 
7 = 2. 

Less trivial are the examples where 

[[H, M],M] = j 2 H + terms commuting with M. (3) 

Note that (jTj) implies Indeed, one can always write [[H,M],M] = >y 2 H + X. 
Using ([T]) there follows 

[X,M] = [[[H,M],M]- 1 2 H,M] 

= [[[F,M],M],M]-7 2 [tf,M] 

= 0. (4) 
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Consider the harmonic oscillator. The Hamiltonian is 

H = huv{btf + 6+6) (5) 

with [b, b*] = 1. Let M be the Hamiltonian of the shifted oscillator 

C 2 

M = H -C{b + tf) + j-. (6) 

Then one finds 

[[H,M],M] = (hw) 2 (H-M) -C 2 hw. (7) 

This is of the form Q with 7 = huj. Note that the choice of counting operator 
for this example differs from the obvious one, this is, M is not equal to b^b. The 
counting operator b^b commutes with the Hamiltonian H and is therefore not 
suited for the present analysis. In the one-parameter family © of allowed counting 
operators only the choice Q = is exceptional because then M commutes with H. 

Also the Jaynes-Cummings model (see for instance |13|. I14j) satisfies ([3]). The 
Hamiltonian is H = Hq + V with 

H = hw{btf + tfb) + ^tw Q a z (8) 
V = hn{b ] a- + ba + ). (9) 
Take M = Hq. A short calculation gives 

[[H,M],M] = h 2 (io -uj) 2 (H - M). (10) 
Hence, (pQ) is satisfied with 7 = K(ujq — uj) (assuming the off-resonance condition 

U) ^ UJq). 

In 1981 Witten [15J introduced his non-relativistic model of supersymmetry. 
The one-particle Hamiltonian reads 

d 2 1 o. , 1 dw 
H = -^d^ + 2 mW ^ + 2 ha ^- (U) 



It can be written as 



with 



where 



H = {Q,Ql} (12) 



Q = — ;= (P — imw(x)) o" + , (13) 
V2m 



h d 1 

idx a± = 2 %(Jy 
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It is straightforward to observe that this supersymmetric Hamiltonian satisfies ([3J) 
with M = o x and 7 = 2. 

A highly non-trivial example of (pQ) is the transverse Ising model [16], with 
Hamiltonian 

h = 4 ? jy<7? °* ( i5 ) 

ij k 

Assume that Jij = Jji and Ja = 0. Choose M = a^. Then one calculates 

[H,M] = 2/^./,,<^. (16) 

[[H,M],M] = 8(H + hM)+4^2j ij a y l a y j . (17) 

u 

The relation between [[H, M],M] and H is not so easy to analyse as in the previous 
examples. However, some further calculation shows that ([I]) is satisfied with 7 = 4. 
Finally consider the Hubbard model. The Hamiltonian is 

H = -J2 t ij Yl h \,° h h° + a Yl n kA n Kl (18) 

i,j °-=t,l k 

with frtj and the creation operators for an electron with spin up respectively 

spin down located at site i, and with nj i(T = b\ a bi >(T the particle counting operator 
at site i and spin a. The coefficients and a are assumed to be real. They must 
satisfy t{j = tji to make H Hermitean. In addition, one assumes that ta = 0. Now 
let 

M = Y. nk ^ nk 'i- ( 19 ) 

k 

Then, calculations similar to those of the transverse Ising model show that (pQ) is 
satisfied with 7 = 1. In this case M counts the number of electron pairs sharing a 
lattice site. 

The above examples make it clear that it is worthwhile to investigate the 
relation ([1]). Its first consequence is that one can write the Hamiltonian H into 
the form 

H = H Te{ + R + R\ (20) 

with H rc { and R satisfying 

[H Ie{ ,M] = (21) 
[R,M] = jR. (22) 
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Indeed, let 



H iei = H — \[[H, M],M], (23) 

7 

R = ^[[H,M],M} + ±-[H,M]. (24) 

It is then straightforward to check that (|20l [2"Tl [2"2"j) are verified. Conversely, if 
a Hamiltonian H can be written as (I20|) , with H re i and R satisfying <\20\ \21\ 1221) 
for some M and 7 then (pTJ) follows automatically. Hence, ([T]) and (f2"Ul [2H l2"2"|) 
are equivalent. The algebraic relation (I22h is also satisfied in Fock space by the 
annihilation operators and the particle counting operator. This is the motivation 
to call M the counting operator. 

If M = E is an orthogonal projection operator then one has 

Hrei = EHE + (I- E)H(I- E), 

R = (I-E)HE. (25) 

This means that H re f is the diagonal part in the block matrix representation de- 
termined by E, while R and R< are the off-diagonal contributions. Note that in 
this case RR) + R^R = -[H, E] 2 . Hence, if [H, E] is a multiple of one then R and 
-Rt satisfy anti-commutation relations. 

A simple calculation shows that H vc i = and R = jhuio(a z + ia y ) in the 
example of the Pauli spin, H re i = M + Q 2 /huj and R = Qb—C^/Tvuj for the harmonic 
oscillator, and H ie i = H$ and R = hnb^a- in the Jaynes-Cummings model. In the 
supersymmetric model is 

R = i^K+i<r„). (27) 
In the transverse Ising model is 

Hr C i = —E-^^ + ^-^E^ < 28 ) 

i,j k 
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R = ^E^H-H-^l) + |E^M- 



(29) 



Note that H TC i is the Hamiltonian of the quantum XY-model. In the Hubbard 
model is 

Hrei = - E ^hhti 1 ~ ( n H ~ n jif) 

i,3 
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k 

R = -^2t ij bl 1 b j ^(l-n il )n jl 

hi 

Another consequence of the assumption (pQ) is that the spectrum of M neces- 
sarily consists of equally spaced energy levels. Indeed, from Mip = flip follows 

M(Rip)=R(M-^)ip = (fi-^)(Rip). (32) 

Hence, either Rip = or /x — 7 is also an eigenvalue of M. Similarly, one concludes 
that either R*ip = or \i + 7 is also an eigenvalue of M. 

Creation and annihilation operators may be used to generate a series of eigen- 
states starting from a single eigenstate, which has been found by other means. 
Let us show that the operators R and R) introduced above can serve the same 
purpose. 

Introduce the notion of a stable eigenstate [TZ] as an eigenstate ifi of H for 
which either Rip = R^ip = or for which numbers x, y exist so that 

(xR + yR))ip = tp. (33) 

The exceptional case that Rip = RJip = is added for convenience. Both eigen- 
states of the Pauli spin Hamiltonian are stable, with x = y = ±2/tjioo. The ground 
state |0) of the harmonic oscillator is stable. Indeed, one has R\0) = — (C 2 /hio)\0). 
Hence, (f33j) is satisfied with x = —ftw/C 2 and y = 0. 

Given a stable eigenstate ip of H with eigenvalue v one can try to construct a 
new eigenstate \ defined by \ = f(M)ip. One has 

H X = [H,f(M)]iP + u X 

= [R + R\f{M)]iP + v X 
= [f(M + 7 )-f(M)]RiP 

+[f(M- 1 )-f(M)]R^ + v X . (34) 

Now let f(u) be a function such that for some real A 

f(M + 7)i?V + /(M - 7 ).Rty 
= (1 + Xx)f(M)RiP 

+(1 + \y)f{M)R^. (35) 

Then the expression becomes, using (|33|) . 

H X = \f(M)(xR + yR^iP + v X 
6 



(30) 
(31) 



(36) 



One concludes that x is an eigenstate of H with eigenvalue u+X. The construction 
of new eigenstates of H is therefore reduced to finding functions f(u) solving the 
eigenvalue equations (|35|) . 

As a first application let us assume that ijj is a stable eigenstate with xy ^ 
and x + y ^ 0. Then there exists at least one other stable eigenstate x °f the form 
X = f(M)ip. Indeed, make the choice f(u) = exp(zit) with z a complex constant. 
Then a sufficient condition for (|35p to hold is 

e zM e z ^Rif; = (l + Ax)e 2M i?V, 
eZ M e - z7 ^ = (x + \ y ) e *M R ^ ( 37) 

These equations have a trivial solution when 

A = -£±» (38) 

The solution is obtained for z satisfying exp(^7) = —x/y. The resulting eigenstate 
X is orthogonal to ip because the assumption x + y ^ implies A 7^ so that ip and 
X have different eigenvalues. It is straightforward to verify that x is again a stable 
eigenstate. However, repeating the above argument starting with x reproduces the 
eigenstate tp. Hence, with this choice of f(u) only one additional eigenstate can 
be obtained. 

When applied to the Pauli spin example the above reasoning allows to derive 
one of the two eigenstates from the other. The value of A is then ±tkVQ. Also all 
eigenvectors of the Jaynes-Cummings model are stable. It is well-known that the 
spectrum of this model consists of a singlet, which is the ground state level, and an 
infinity of doublets. The ground state is stable in a trivial way because it satisfies 
Rip = R'lp = 0. Each of the doublets consists of two stable eigenvectors which 
transform into each other by the above mechanism. In the one-dimensional Hub- 
bard model with 4 sites and a half filled band one can show [S] that several pairs of 
stable eigenstates occur which transform into each other by the above mechanism. 
In addition, these states are eigenstates of the anti-commutator {R, R^}. Similarly, 
in the one-dimensional transverse Ising model with 3 and that with 4 sites one can 
show that all eigenstates are stable and that they either satisfy Rip = R^ip = or 
they occur in pairs which transform into each other in the way described above 

IB]- 

The ground state of the harmonic oscillator does not satisfy the condition 
xy 7^ 0. Hence, the previous result cannot be used to construct a second stable 
eigenstate, in agreement with the observation that the ground state is the only 
stable eigenvector of the model. In this case, no new eigenstates of the form 
/(M)|0) can be obtained from the ground state |0), as can be seen immediately 
using the commutation relation [M, b] = —hwb. However, it is well known that all 
eigenstates are obtained by repeated action of the creation operator &t = Rs /(, + 
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More generally, if in some model an eigenvector ip of H exists such that xRip = 
ip then an annihilation operator B, annihilating ip, can be defined by B = R — x~ l . 
However, in general not much is known about the commutators [B, B^] = [R, R^] 
and [H Ief ,B^]. For this reason, no further progress was made in this case. 

Once a pair of stable eigenvectors has been obtained one can try to find func- 
tions /, other than exponential ones, satisfying (|35[) , This does indeed work for 
the following artificial example. Fix k, £, /i, and v so that k / 0, |i / 5) an d 
£ + v = 2/j>. Let 

H = I k fi K I (39) 



and 





M = . (40) 



The eigenvectors and eigenvalues can be calculated explicitly. Two of the eigen- 
vectors are stable. They are of the form ip = (a 2 ,2a,2) T . The third eigenvector 
is of the form \ = (1><5j — 1) T and satisfies x = f{M)ip, with /(«) of the form 
f{u) = A + Bu. One concludes that the three eigenvectors together form a triplet. 
Note that this example of 3-by-3 matrices has been considered in [11] as well. One 
can write H = H TC { + kV and I = M + kW, with kV = R + R^ . However, there 
does not exist a ^-independent matrix W such that H and / commute for all values 
of k. Hence, in the terminology of [TT], the pair (H,I) is not integrable. Similar 
triplets exist in the 4-site Hubbard ring — see [6] . 

In summary, we have studied model Hamiltonians H together with a Hermitean 
operator M such that the third commutator of H with M is proportional to the 
first commutator. Then the spectrum of M consists of equidistant levels and 
M is called the counting operator. In this context one can define operators R 
and i?t which have some similarity with annihilation and creation operators. An 
eigenvector ip of H is said to be stable if either Rip = R^ip = or there exist x, y 
such that (xR + yR^)tp = ip. If xy ^ and x + y ^ then a new stable eigenvector 
can be constructed whose eigenvalue is decreased with (x + y)jxy. This is the 
main result of the present paper. But the analysis also suggests to decompose 
the spectrum of H into multiplets determined by the symmetries of the reference 
Hamiltonian H rci = H - R- RK This view is supported by the detailed analysis 
of some small systems, including the 4-site Hubbard ring, which is studied in [B]. 
Acknowledgements We are grateful to Prof. J. Perk for pointing out [T91I201 
and to Prof. M.S. Plyushchay for pointing out |22j . 
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